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ABSTRACT 

The reception of FM carriers is characterized by a capture effect whereby the 
message of the dominant carrier is recovered when two or more FM carriers are 
present. This research uses a computer to form a type of average of the instan- 
taneous frequency of the receiver input. The results establish that averaging 
(smoothing or lowpass filtering) of the instantaneous frequency reveals the capture 
effect. The effect on capture of bandwidth and amplitude ratio for the case of two 
carriers is revealed. The results show that capture can occur when two Carriers are 


separated by as little as 0.17 dB in amplitude. 
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I. INTRODUCTION 


Radio communication usage continues to expand because the military and 
commercial markets seek (1) personalized service such as paging and cellular 
applications, (2) ever greater exchanges of text and data (facsimile, computer 
networking), and (3) transfer of high quality signals (television). The messages are 
"traffic". The carriers are sine waves. 

There are two primary ways of carrying messages using sinusoids. The message 
can cause to vary (modulate) the amplitude or the argument (angle) of the sinusoid. 
These are called amplitude modulation (AM) and angle modulation (4M) 
respectively. 

Amplitude modulation can assume various forms such as double sideband with 
carrier, double sideband suppressed carrier and single sideband. Angle modulation 
can be either frequency modulation (FM) or phase modulation (PM). With analog 
messages, there is little difference between FM and PM. In fact, FM hardware 
would be used should PM be chosen. Consequently, PM with analog messages is, in 
practice, reduced to FM. [ Ref. 1: pp. 198-199 | 

Frequency modulation exhibits several interesting characteristics not shared by 
AM. Two such characteristics are threshold effect and capture effect. The threshold 
effect relates to the quality (signal—to—noise ratio) of the demodulator output as a 
function of the quality of the received signal (demodulator input). Below 
"threshold" quality of the input, the output quality deteriorates rapidly. The 
capture effect relates to the ability of the demodulator to recover the message of the 


dominant carrier when two or more FM carriers are present. 


This research is concerned with the capture effect of FM. Computer simulation 
is used to verify an idea concerning the reason for (origin of) the capture effect of 
FM demodulators. 


Background information on FM and the capture effect is presented in chapter Il 





of this thesis. Chapter III describes the computer simulation procedure and 
program. Results of the simulation are contained in Chapter IV. Conclusions and 
recommendations are listed in chapter V. The computer program is listed in the 


Appendix. A list of references is provided. 


Il. BACKGROUND 


A. EQUATION FOR FM CARRIER 


An angle modulated sinusoid is expressed as 


s(t) = B cos|¢,(t)] (221) 
where the constant B is the carrier peak amplitude. A complete oscillation cycle 
occurs whenever the argument ¢,(¢) changes by 27 radians. The instantaneous 
frequency of the angle modulated wave s(t) is defined by 

MO) = oat aa 
Frequency modulation is that form of angle modulation for which the 
instantaneous frequency ff t) varies in accordance with the baseband signal or 


message voltage m(t). Typically, a linear relationship is chosen so that 


(AU ae k em{t) (2.3) 
where f is the frequency of the unmodulated carrier and the constant / f represents 
the frequency sensitivity (in Hertz per volt) of the modulator. From Eqs. (2.2) and 
(2.3), the angle d {t) is 


b(t) = 2nft+ dak pf 'm(r) dr (2.4) 


where, for convenience, the angle of the unmodulated carrier wave is assumed to be 


0 at ¢ = 0. The frequency modulated wave is, therefore, 


nn 
to 
Cyr 

ee” 


s(t) = B cos[27-f,t + dak pf 'm(r) dr 


B. CAPTURE EFFECT 
Consider the case of two FM carriers, s;(¢) and so(¢), present at the demodulator 
input and in the same frequency band. It is well known in practice that the output 


of the demodulator is the message of the dominant carrier. That is, if 


si(t) = By cos|[2nf,t + dak pf m(r)ar 
so(t) = Bo cos[2a( f+ e)t + dak f °ma(7) dr (2.6) 


then the output of the demodulator is proportional to 
m(t), when By> Bo 
ma(t), when B;< Bo 


The carrier frequency offset € is any value such that both s;(t) and so(t) are in 
the same operating frequency band of the receiver/demodulator. This inherent 
ability of an FM demodulator to suppress weaker FM carriers is called capture 
effect. 

The FM demodulator can be a phase locked loop, a slope detector which 


includes Foster—Seeley and ratio detectors, or a pulse counting discriminator. 


C. ANALYSIS 

As shown by Eq. (2.3), the instantaneous frequency ff t) of a frequency 
modulated carrier is proportional to the message m(t) being carried. Consequently, 
it is claimed in the literature that the output voltage of a frequency demodulator is 
proportional to f(t) — f= Af{t) [ Ref. 2: pp. 169—172 |. This is indeed the case in 
practice when the demodulator input is a single frequency modulated carrier. 


However, it is not the case, because of the capture effect, when the demodulator 


input is multiple frequency modulated carriers. 
The instantaneous frequency of the sum of two frequency modulated carriers 


can be derived as follows. Let 
v(t) = By cos¢;(t) + Bo cosdo(t) (2.7) 


Using trigonometric identities, Eq. (2.7) can be written as 


ut) = Ccosd¢s(t) (2.8) 
where, 
CAG) = By2+ Bo2 + 2B, Bocos| ¢;( t)—¢2(t)] (2.9) 
Bysing,(t) + Bosingoe(t) 
mi) tists (2.10) 


Bycos¢,(t) + Bocos do(t) 


From Eq. (2.2), the instantaneous frequency f(t) is 


glue d tan”! Bysin 1 t)+ Bosin 2 i] 
Qn dt 1cOS ¢;(t)+ Bocos do{t 
_ (Bi ¢1'cos d;+ Bod2'c0s b2 ) C+( Bi ds'sin G1 + Bo Go'sin G2) D 


2m(C2+ D2 ) 


_ Bi2d;'+ Bo? G2'+ | Bi Bo c08 ( b1-G2 )| ( G1'—69' (2.11) 
2 7| By2+ Bo 2+2B, Bocos( ¢1-¢2)| 


where, 
C= B, cos¢;(t)+ Bo cos¢o(t) 
1B = By sin go(t)+ Bo sin b(t) 


and ¢;', ¢9' represent £di(1) and £d2(t) respectively. 


As before, for an FM carrier 
b(t) = 2aft+ Ink, f'mA{r)dr 
t G i 9 ? 


b(t) = 2nf, + 2akem{?) i= 1,2 I) 
where the carrier frequencies for both messages are assumed to be the same and 


mf t) = 1 th message 


k > frequency sensitivity of the frequency modulator in Hertz/volt 


Equations (2.8) through (2.12) can be reduced and interpreted for particular 


messages. For example, let the two messages be tones where 


m {t) = ID cos27 f(t) alee (2.13) 


For simplicity, let D, = Dy = 1. Then, 
b(t) = 2nftt 2k, f 'cos(2af-7)dr 
1 C f 0 t 
k 


- anit =p sin2a fit 
2 


g' {t) = nf, omk ecosanfit lee (2.14) 
Using these values in Eq.(2.11) gives 


FAB 24 B2)4 ‘X+B,BoY:-Z 


OTS (25) 
B,2+Bo2+2B, BoZ 


where 
X = By? cos27ft + Bo2 cos2zrft 
e— 2f,+k -(cos2aft + cos2z ft) 


= : ae ee 
ZL = cos|2mh e( aT sin2afot + Dail sin27f t )| 


Equation (2.15) is plotted in Figs. 2.1.a through 2.2.c along with m;(t) and 
mo(t) where f; = 1 Hz and fo = 1.3 Hz. It can be seen that there are many 
impulselike frequency deviations from the message waveforms. The variable in Fig. 
2.1 is 6 ( modulation index for FM). The bandwidth of the modulated carrier is 
related to f. In general, bandwith increases linearly with 6 [ Ref. 3: pp. 194-197]. 

There are two major factors that affect the graph of instantaneous frequency. 
First, the bandwidth of the modulated carrier affects the spacing in time between 
adjacent impulse—like frequency deviations. Figures 2.1.a through 2.1.c show that as 
the bandwidth decreases, the spacing increases. Second, the difference in magnitude 
of the two modulated carriers affects the extent of the deviation. Figures 2.2.a 
through 2.2.c show that as the magnitude difference increases, the deviation extent 
decreases. 

However, these results do not show specific evidence of capture. Rather, they 
suggest that the instantaneous frequency needs to be averaged over a short time 
interval to reveal the capture effect. The same thing is true for the case of close (but 
different) carrier frequencies, as shown in Chapter III. 

The idea that the capture effect is caused by the lowpass filtering operation of 
frequency demodulation led to this research. Since lowpass filtering can be thought 
of as averaging (or smoothing), it was decided to investigate the capture effect by 
averaging the instantaneous frequency using a computer. 

So, in hardware, the averaging is done by the lowpass filter (LPF). A LPF is 
part of every frequency demodulator (frequency—to—voltage or voltage—to-frequency 
converter). Consequently, the output of a frequency demodulator is a lowpass 


filtered version of instantaneous frequency deviation A ff ie 


Again, lowpass filtering can be considered as a smoothing operation on the 
applied voltage. In this research, we used computer simulation to investigate the 
effect of smoothing of A fi t) to reveal the presence of capture. Smoothing is created 
by averaging time intervals between adjacent zero crossings of the FM carrier. Since 
this interval (7) represents the time required for 7 radians excursion of the 
argument of the cosine carrier, and since 7 radians is taken as a half cycle, then two 
such intervals (7; + Tz) become a full cycle and Tet, is the number of cycles per 
second or the instantaneous frequency for 7; + Jo seconds. Averaging N such 
groupings provides a smoothed value of or short-term average of the instantaneous 
frequency over that expanded time interval. 

A computer program was used to approximate 7; and 7> for every interval 
between zeros of v(t). These values are then averaged and the results are plotted. 


The program is explained in Chapter III. 
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HI. SIMULATION PROCEDURE 
The FM demodulator, including LPF, is simulated by a computer program so 
that critical parameters, such as signal amplitudes, bandwidth of the FM system, 
message forms, and interval of averaging, can be easily changed. Since the sampled 
version of the demodulator input forms a sequence of numbers, it is convenient to 
consider it as a row (or column) vector. Once that is done, the rest is simple 


manipulation of the row (or column) vector. 


A. WAVEFORM GENERATION 


The desired waveform v(t) can be generated using a simple routine 
(WAVEGEN). In this routine, the carrier amplitudes, values of k p carrier frequ- 
encies and message frequencies are input parameters which can be changed. This 
resulting wave is sampled at a rate which will provide the desired accuracy in 
measuring the time between adjacent zeros of v(t). In this work, a carrier frequency 
of 200 Hz and sampling rate of 12 kHz are used. 

The sampled values contain no information about frequency. Rather, the 
number of samples between zero crossings are related to frequency because the 
number of samples between zero crossings itself represents a half cycle oscillation of 
argument of v(t) and its reciprocal can be interpreted as twice the frequency of that 
period. Therefore, the routine notes the polarity of each sample so that the number 
of samples between zero crossings (polarity changes) can be easily counted. This 


result is then used as the input for the next routine. 
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B. COUNTING 

The number of samples between zero crossings (polarity changes) multiplied by 
the fixed sampling interval (A) is considered to be taken as half the instantaneous 
period of that interval. So, a sequence of integers which consists of the number of 
samples between zero crossings need only be generated. The 7 th entry for this 
sequence is denoted by Yr. , t= 1, 2,3,---+ as shown in Fig. 3.1. The totality of 


entries forms a sequence of integers. 


u(t) / voll 


ip Ve ule t/sec 


EEE ER EERE EEE EEE EEE EY EEE EEE RE E> 
AT AT3 


Fig. 3.1 Formation of sequence 


Furthermore, r itself is a time index. Since there are 2— 1 entries before the 


current interval r,, T, is the period for the time iD, liye sf, . (The middle of 
the interval r, is selected for the time index of interval r, . For example, the time 
index of T3 is 7) + T2 + 573.) 

Since the first and last interval may contain a number of samples that does not 
correspond to a complete interval between adjacent zeros, these two entries are not 


used in the averaging routine. 
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C. AVERAGING 


Given the sequence of integers, any average T), over N intervals is formed as 
k+N 


k = any positive integer (3.1) 





r, = ? th interval 


N = number of intervals averaged 


Since frequency is taken as the reciprocal of the period, the averaged frequency, 





denoted by Jj, , is given by 


= (3.2) 
2T,, A 





I 


To obtain more realistic (smoother) results, a moving average is used rather 
than fixed time interval averaging. To illustrate this, assume a sequence of intervals 


as follows 
imi a ea 


Then, if the average is done over 3 intervals, the averaged version of the sequence 


: will have the following entries 


ie = T+ “34 T3 
T» = To+ sgt 
T+ 7 + T,, 
T= fe ora 5 (3.3) 
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Again from the relationship between period and frequency (Eq.(3.2)), averaged 
version of frequency will form the sequence 


i ] ] 


2T A 2T2r4 = 2T,r 





In this case, the time indices for the averaged version of intervals are 


ee = time index for T; 


oe = time index for 7» 


T+ gy = time index for T3, and so on. 


Generally, time indices for the moving average over N intervals are found as 
i—] ) +1 >. : 
pipet pu; Ty (time index for T, ) (3.4) 
Since the values of each 7} are not the same, the averaged frequency values when 


plotted will not be equally spaced. 
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IV. SIMULATION RESULTS 
With an average of 2 intervals (N = 2), the resulting plot is very similar to that 
of the analytical one. That is, when almost no averaging is done, the resulting plot 
should yield that of the analytical one. This is another way of verifying that the 
simulation program is running properly. As the number N of average intervals 
increases, the impulse—like frequency deviation becomes smaller and the overall 
frequency deviation approaches that of the stronger signal. 
Figures 4.1.a through 4.5.c are the simulation results. The results are presented 
as averaged instantaneous frequency versus time. In all cases, the message m(t) 
corresponds to that of the dominant carrier. The following variables are used to see 
the effect of system parameters. 
a,/ a2 = ratio of carrier amplitudes (1 through 1.5) 
6 = modulation index (1 for narrowband, 10 for wideband, 50 
for very wideband FM) 
fi = message 1 frequency (1 Hz for sinusoidal message) 
fp = message 2 frequency (1.3 Hz for sinusoidal message) 
f= carrier frequency (180 Hz through 200 Hz) 
N = number of intervals averaged (2, 20, 40, 60, 80) 


Various messages were used in the simulation. Figures 4.1.a through 4.1.c show 
the results when the messages are linear intrapulse FM (up CHIRP and down 
CHIRP). Figures 4.2.a and 4.2.b are for exponentially decaying and growing 
sinusoidal messages. Sinusoidal messages are used extensively to show the results for 


various correlations of the listed variables. 


Li 


Even though there may be other factors which influence the capture effect, only 


the following four factors are investigated in this research. 


A. EFFECT OF N (AVERAGING TIME) 

Figures 4.3.a through 4.3.d shows the effect of N. The other factors (a1/a2, 3, A, 
fh, fe) remain fixed. This factor is essentially the cause of the capture effect. As N 
increases, the instantaneous frequency increasingly resembles the message of the 


stronger signal. 


B. EFFECT OF SIGNAL AMPLITUDE RATIO 

Figures 4.4.a through 4.4.d show the effect of signal amplitude difference. While 
the other factors remain fixed, a;/a2 varies from 1.01 to 2.0. The effect of signal 
amplitude difference is very interesting. Even with very small signal amplitude 
difference, the system still locks onto the stronger signal. For numerical comparison, 
the 0.8 dB difference case (Fig. 4.4.b) and 1.6 dB case (Fig. 4.4.c) show little 
difference. In the case of very small amplitude difference, more averaging is needed 
to see the capture effect. Figure 4.4.a indicates capture for a carrier separation of 


0.17 dB (ay/GpeeleO2) 


C. EFFECT OF MODULATION INDEX 

Modulation indicies are changed from 1 to 50 in Figs. 4.5.a through 4.5.c. 
Capture occurs for both narrowband FM and wideband FM. Since the narrower 
band FM shows fewer impulse—like frequency deviations, each of longer duration 
than those of wider band FM, narrow band FM needs to be averaged longer (see 


Fig 23e 
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D. EFFECT OF CARRIER FREQUENCY DIFFERENCE 

Figures 4.6.a and 4.6.b show the effect of carrier frequency difference on 
capture. It appears that capture becomes more complete as the carrier frequencies 
separate. Obviously, the two carriers must lie in the passband of the receiver 
(intermediate frequency amplifier) if they are to be present at the input to the 


demodulator. 
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where a@;/ a2 = 1.5 
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V. CONCLUSIONS AND RECOMMENDATIONS 

The simulation results establish that the lowpass filtering portion of frequency 
demodulation accounts for the capture effect of FM receivers. Without filtering, the 
instantaneous frequency deviation of the sum of two carriers has an impulse—like 
history which does not consistently favor the message of the dominant carrier. With 
averaging or smoothing (lowpass filtering) of the instantaneous frequency deviation, 
capture occurs. 

The simulation results indicate capture occurs for both narrowband and 
wideband FM. Capture is evident when the carriers are separated in amplitude by 
less than 0.2 dB. The effect of averaging time is revealed by the simulation output. 
In general, longer averaging times produce smoother approximation to the message 
of the dominant Carrier, as expected, and also show the capture effect more 
convincingly. 

It is recommended that the capture effect be verified using an operating 
experimental system in which system parameters can be controlled and accurately 


measured. 
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APPENDIX 
Program List 

This appendix contains the program used to calculate f, defined in the thesis. 
The program actually determines T. and T,. as given by Eq. (3.1). 

The input equation is s;(t) + so(t) as defined by Eq.2.6. This equation is 
sampled 12,000 times in the one second simulation interval. The interval between 
samples is, then, 1/12,000 = 83.3 usec = A. The polarity of each sample value is 
noted. The number of consecutive sample values of like polarity in interval 2 is r, 

The program stores each value of 7, and forms T,, from Eq. (3.1) for various 
assigned values of N. Then h is calculated according to Eq. (3.2). Finally, h. is 
plotted for each time index of Eq. (3.4). 

The program contents are following: 

MAIN : main routine for the simulation of capture effect 

CHIRPGEN _ : subfunction routine for generation of CHIRPs 

WAVEGEN - : subfunction routine for generation of sinusoidal waves 

COUNT : subfunction routine for counting and averaging 

SAMPLE : subfunction routine for display 
The program was written in APL and run on IBM PC/AT with 3 Mb of memory. 
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MAIN 


V Y¢MAIN A3SA;Y; TNFPSMP SESE Ee 


own routines. 
The output will be, 


(not equally spaced) 


frequency 


D> 2 D> D> > 2 2 22 2 2 


O«+'Generating waves....' 


A Generating arbitrally sinusoidal waves 


TMPSMP¢WAVEGEN A 
U0 Count tiger en 


nA Count the number samples between 
A zero crossings and form a resulting sequency 


YECOUNT TMPSMP 
Ue" Samp linge 


This 1S main routine for the capture effect 
Simulation. For the clality of logical flow, 
each simulation step is divided into their 


T: time index fopreelrancuen 
TIME: time index for aceraged frequency 
F1,F2: individual instantaneous 


RES: averaged version of frequency 


\ 


A Since the numbers of elements of sampled time and 
A message are too many for plotting purpose, it is 
A sampled without missing any important information. 


T¢SAMPLE T 
FI¢SAMPLE Fl 
F2¢SAMPLE F2 
licGabone. ' 
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CHIRPGEN 


Vv 


De D D Dep 


VOCCHERPGEN AM MIO; ALSA2; PHI; PH2;V ;iMPL Ime 


This routine generates two CHIRP signals, namely up 


and down CHIRP to test the program. 
Equations for this routine are, 
F1=100T+150 for up CHIRP 
F2=250-100T for down CHIRP 
where 0<T<1 
OTO<0 
A2¢1 


TSLICE¢i (€1.3%0.5)x10000 
met tliokler )=TSLICE 
D€50 
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TMP¢TMPXT 
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TMP¢TMPXxT 
TMP1€(02)xD 
TMP1¢TMP1XT 
TMP1¢TMP1XT 
PH2¢TMP-TMP1 
TMP¢20PH1 
TMP1¢20PH2 — 
TMP¢A1XTMP 
TMP1¢A2XTMP1 
V¢TMP+TMP1 

VexVv 

F1€¢2xDxT 
F1¢F1+200-D 
F2€¢200+D 

TMP¢2xD 
TMP¢TMPXT 
F2¢F2-TMP 
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A 


A Routine to generate samples of sinusoid waveforms 


A 
DIO<0 


A Second wave amplitude. This routine accepts the first 
A signal amplitude as an input argument 


A2¢1 
P2€902 


TSLICE¢l (1.540.5)x10000 
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TMP1¢P2xXFC1xT 
TMP¢TMP+TMP1 
V1¢20TMP 
V1¢V1xAl1 
TMP¢P2xFM2 
TMP¢+TMPXT 
TMP¢+TMP+PH2 
TMP¢10OTMP 
TMP<¢TMP-(1°0PH2) 
TMP¢TMPXKF2+FM2 
TMP1¢P2xFC2 
TMP1¢TMP1xT 
TMP¢TMP+TMP1 
V2€¢20TMP 
V2¢A2xV2 
Vex(V1l+v2) 
TMP¢P2xFM1 
TMP¢TMPXT 
TMP¢TMP+PH1 
TMP¢KF1XCOSCTMP) 
F1¢FC1+TMP 
TMP¢P2xFM2 
TMP¢+TMPXT 
TMP¢+TMP+PH2 
TMP¢€KF2xXCOS(CTMP) 
F2¢FC2+TMP 
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D. 


C1] 

f2) 

C3] 

C4) 

C5] 

C6] 

C7) 

C8] 

C9] 

C10] 
C11] 
C12] 
C13] 
C14] 
C15] 
64 
tae J 
C18] 
bod 
C20] 
2] 
b22) 
h23) 
C24] 
C25] 
C26] 
C27] 
C28] 
C29] 
C30] 
ees 
C32] 
C33] 
C34] 
C35] 
C36] 
C37] 
C38] 
[39] 
C40] 
C41] 
C42] 
C43] 
C44] 
C45] 
[46] 
C47] 
C48] 
C49] 
C50] 
C51] 
C52] 
(53) 
[54] 


COUNT 


V Y€COUNT V;C;COUNTER;s INIT; I; LIMIT;:D;:Y;TMP 


; EXTIME1 ; EXTIME2 ; INTV 


ARoutine to form half period sequence and average it 


A 
A Set the initial matricies and values 
A 
Ceg 
COUNTER¢€0 
INIT¢€VCOIO] 
I¢0IO 
LIMIT¢(pV)+01I0-1 
A 
A Routine to count number of samples between 
A Zero-crossings 
fe) 
LOOP1: 
*C(VLCIJ=INIT) /LOOP2 
C¢+C, COUNTER 
INIT¢V(TI] 
COUNTER€1 
HOOr 2. 
LOOP2: 
COUNTER¢COUNTER#+#1 
LOOP3: 
I¢I+1 
+(I=LIMIT)/LOOP1 
Cee, COUNTER 
A 


A Drop the incomplete cycle at the beginning and the end 


\ 


A 

Bor iIMereclOLOI<(1l=TSLICE) 

Cce1!lc 9 Ce711lCc 

D¢e§ 

TIME¢9 

fs) 
A Average with the value of INTV 
A 

INTV¢60 

EH LO 

LIMIT¢€(eC)-COUIO+1+INTV) 

A 

A Calculate extra time due to averaging 
A 

Pee tME2eC+/INTVtC)+TSLICE 

e 

A Calculate the time index which is not equally placed 
A 

LOOP4: 

TMP¢+/INTVTC 

D¢D,TMP+INTV 

TIME¢TIME,CCOIO] 

Cearec 

eit 1 
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[55] seb MIT )7 EOOr4 
[56] TIME¢((+\TIME)+TSLICE) tEXTIME1+ CEXTIME2+2 ) 
C57] A 
[58] a Calculate averaged frequency for the time indicies 
C59] A 
C60J Y€¢C1+2xD)xTSLICE 
Vv 
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E. SAMPLE 


Veet ©SAMP LEeexs 0 TOs; SdZejeloTiMelT ; dey X 


Ci) 

[2] a Routine to sample for plotting (number of samples = 60 
C3] 

C4) OIO¢o 

tS) Y¢® 


[6] SIZ¢l (CeX)+600) 

(7] LIMIT¢€(CpX)+0I0-1 

C8] I¢0IO 

[9] LABEL1: 

C10] Y¢Y ,XCIXSIZ) 

[11] T¢eI+1 

C12) 2C€CIXSIZ)<LIMIT)/LABEL1 
Vv 


3f 
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